Keck Observations of the Galactic Center Source G2: Gas Cloud or Star? 
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We present new observations and analysis of G2 - the intriguing red emission-line 

(SJ ', object which is quickly approaching the Galaxy's central black hole. The observations 

were obtained with the laser guide star adaptive optics systems on the W. M. Keck I 

<T^ and II telescopes and include spectroscopy (R ~ 3600) centered on the Hydrogen Br- 

7 line as well as K' (2.1 fim) and V (3.8 ^m) imaging. Analysis of these observations 

shows the Br-7 line emission has a positional offset from the V continuum. This offset 

♦pH, is likely due to background source confusion at L' . We therefore present the first orbital 

solution derived from Br-7 nne astrometry, which when coupled with radial velocity 

measurements, results in a later time of closest approach (2014.21 ± 0.14), closer pe- 

riastron (130 AU, 1900i? s ), and higher eccentricity (0.9814 ± 0.0060) compared to a 

solution using L' astrometry. The new orbit casts doubt on previous associations of G2 

and a low surface brightness "tail." It is shown that G2 has no K' counterpart down to 



(N 
> 



K ~20 mag. G2 s L continuum and the Br-7 line-emission is unresolved in almost all 
epochs; however it is marginally extended in our highest quality Br-7 data set from 2006 
and exhibits a clear velocity gradient at that time. While the observations altogether 
suggest that G2 has a gaseous component which is tidally interacting with the central 
black hole, there is likely a central star providing the self-gravity necessary to sustain 
the compact nature of this object. 
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1. Introduction 



Recently, iGillessen et al.l (|2012l . 120131 ) reported the discovery of G2, an extremely red object 
with spatially resolved Br-7 emission. The object was interpreted as a ~3 Earth mass gas cloud 
based upon an inferred low dust temperature, observed elongation of the Br-7 emission along the 
object's direction of motion, and a claimed tail along the same orbital trajectory as G2. This 
interpretation is particularly interesting because G2 is on a highly eccentric orbit with closest ap- 
proach to our Galaxy's central black hole within the next year, potentially allowing us to observe 
an unprecedented accretion event on to a supermassive black hole and offe r insight into the re- 
gion surrounding t he black hole (e.g. 



Saitoh et al. 



2012; 



Bartos et al. 



2013 



Morris! [20I2I: iMoscibrodzka et alJboid ; lAnninos et~all bold : 



Yusef-Zadeh k Wardlel[20li ). 



The interpretation of G2 as a gas cloud, however, is not definitive. One challenge for the pure 
gas cloud scenario is that, given the strong tidal fields in this region and G2's low s elf-gravity, G2 
needs to have formed q uite r ecently (~1995, just prior to the initial observations; iBurkert et al. 



2012 



Schartmann et al. 



20121 ). Since such a gas cloud would be tidally disrupted during its periapse 



passage in the upcoming year, the gas cloud model implies G2 will be observed over almost the 
exact extent of its e xistence. Therefore se v eral alternative scenarios invoking an underlying sta r 



have been proposed (JMiralda-Escude 



2012 



Murrav-Clav fc Loeb 



2012 



Scoville fc Burkertll2013T ). 



In these scenarios, the observed heated gas is posited to be circumstellar, either intrinsic or a 
consequence of interaction of the star and surrounding ambient gas. In these stellar scenarios, 
G2 not only has existed for timescales much longer than the observed time baseline, but will also 
survive periapse passage. 

Regardless of its nature, G2's properties and possible origin depend critically on its orbital 
parameters. These parameters have been estimated from observations wi th a short time baseline 
compared to the orbital period (~10 vs.~200 years; Gillcsse n et al.ll2013r) and in a very crowded 



region, making the orbital solution susceptible to biases (JHartkopf et al. 



2001 



. We therefore present 

new measurements and analysis of G2 that minimize the effects of source confusion on estimates 
of G2's orbital parameters and examine the temporal evolution of G2's properties. 



2. Observations 



Two types of new data were collected for this st udy using the laser-guide-star adaptive optics 
(LGS AO) systems at the W.M. Keck Observatory (jWizinowich et al.ll2006l : Ivan Dam et al.l 120061 : 



Chin et al.ll2012f ). Spectra were obtained using the OSIRIS integral field spectrograph (ILarkin et al 
2006) through the narrow-band Kn3 filter, which is centered on the Br-7 hydrogen line (2.1661 /im), 
at a spectral resolution of R ~ 3600. Imaging data were obtained in the i^'-band filter (Ao = 
2.124 fj,m) and the L'-band filter (Ao = 3.776 /im) using the Keck II near-infrared camera (NIRC2, 
PI: K. Matth ews). These data were obtained and reduced in a sim ilar manner as in our previous 



publications (JLu et al. 



2009 



Yelda et al. 



2010 



Meyer et al.1 120121 ) and the specific spectroscopic 



and L' observations utilized in this paper are described in table [Tj 

3. Analysis 

3.1. OSIRIS IFU Measurements 

Because G2 has no detectable K continuum and is fainter than any object we have pre v iously 
extr acted, some ana lysis steps differ from the spectral extraction presented in iGhez et al.l (J2008) 
and iDo et al.l ()2009l ). For all epochs, we created a combined data cube before extracting G2's 
spectrum rather than extracting spectra from individual cubes. The OH sky lines in the data are 
subtracted using sky frames scaled to the strength of families of OH lines in the observed frames 
to account for temporal variations. In 2012, the Br-7 emission line from G2 is coincident with the 
2.180 ^m OH line, so to minimize the systematic effects associated with OH line subtraction, we 
scale the sky only to this line. 

An iterative process was required to estimate the position and spectral properties of G2 in the 
OSIRIS data cubes (see figures [Tj and [2]) since G2's position is needed to place the aperture for 
spectral extraction, and G2's spectral properties are needed to determine which OSIRIS channels 
should be used to measure its position. G2's position was first obtained by visual inspection of 
the data cube. Then, an initial spectrum was extracted at this position using an aperture with 
a radius of 35 mas. Emission from local ambient gas was subtracted using a region free of stellar 
halos within ~0 / /5. A Gaussian was fit to the resulting emission line. In order to refine the position, 
the three-dimensional data cube was median collapsed over the wavelength range corresponding to 
twice the standard deviation of the Gaussian fit to the emission line, centered on the line peak. 
Continuum emission was subtracted by averaging the median of the 25 nearest channels on either 
side of the line. G2's position was further refined with a two-dimensional Gaussian fit to the 
continuum-subtracted image. The iterative extraction process was repeated again to obtain a final 
position and spectrum for G2. Further iterations produced no significant changes. Measuring G2's 



position on a continuum-subtracted frame (see figure 2(b) ) removes the effect of source confusion, 
thereby avoiding astrometric biases. 

G2's spectral properties were obtained using the final Gaussian fits to the Br-7 emission line. 
The radial velocities (RVs) for G2 were calculated from the offset of the line from the rest wavelength 
(A# r7 = 2.1661 yum) and corrected to the local standard of rest. The reported FWHM measurements 
were corrected for instrumental broadening (FWHM=85km s _1 ). Line flux measurements were 



Table 1. Summary of Observations and Measurements of G2 
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AO type/ 
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^'frames 


^ "frames 




FWHM 


Orig 
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Br-7 FHWM 


G2 ARA b 


G2 ADcc 




Date 


Telescope 


(mas) 
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used 




(mas) 


Pub. a 


(kms -1 ) 


(kms -1 ) 


(mas) 
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OSIRIS, 


Kn3 












2006 Jun 18,30; Jul 1 


2006.495 


Keck II LGS 


35 


28 


27 




74 


3 


1125 ±9 


135 ± 30 


222.97 ±5.05 


-105.40 ±2.40 


2008 May 16; Jul 25 


2008.487 


Keck II LGS 


35 


22 


21 




78 


4,0 


e 


e 


181.82 ±6.90 


-64.97 ±2.51 


2009 May 5,6 


2009.344 


Keck II LGS 


35 


21 


19 




79 





1352 ± 32 


171 ±121 


176.26 ± 2.74 


-65.76 ±2.04 


2010 May 5,8 


2010.349 


Keck II LGS 


35 


17 


16 




82 





1479 ± 29 


262 ± 165 


166.29 ±9.11 


-43.95 ±11.05 


2012 Jun 9,11; Aug 11,12 


2012.613 


Keck I LGS 


20 c 


27 


21 d 
NIRC2, 


L' 


68 





2071 ± 146 


706 ± 250 


103.16 ±4.07 


-8.79 ±11.36 


2003 Jun 10 


2003.440 


Keck II NGS 


10 


12 


12 




86 









260.1 ±4.5 


-154.6 ±5.1 


2004 Jul 26 


2004.567 


Keck II LGS 


10 


11 


11 




80 


1 






262.6 ±4.8 


-140.4 ±4.2 


2005 Jul 30 


2005.580 


Keck II LGS 


10 


62 


56 




87 


2 






251.6 ±3.4 


-129.0 ±2.6 


2006 May 21 


2006.385 


Keck II LGS 


10 


19 


19 




81 









226.2 ± 1.9 


-99.4 ±0.8 


2009 Jul 22 


2009.561 


Keck II LGS 


10 


4 


4 




85 









175.4 ±4.9 


-70.7 ±2.1 


2012 Jul 20 - Jul 23 


2012.562 


Keck II LGS 


10 


1316 


1314 




93 









105.9 ±0.8 


-211.1 ± 0.6 
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a References: (0) This Work 


; mlGhez et alj J2005h; (2)1 


Hornstein et al.1 J2007t); 


(3) iGhez et al 


.1 J2008h; 


(4) iDo et aJ 


. (2009) 





b Offset from SgrA*-radio. 

C A smaller square dither pattern of 0'.'5 X tf.'5 was used. 

d In 2012, a more stringent quality cut (FWHM < 68mas) was used because G2 is closer to stars than in other epochs and because the improved LGS performance on Keck I allowed 
this more stringent cut. 

°In 2008, low SNR and changes in the local standard of rest velocity between the two observation dates prevent reliable line measurements. 



made by comparing the Gaussian fit to the non- variable star SO-2 between 2.17-2.18 [im. Statistical 
uncertainties were taken from the root mean square (rms) error of Gaussian fits to three independent 
data subsets; however, we note that systematic errors (e.g. OH sky lines) are not included. 

The final Br-7 line maps were used to evaluate G2's physical size and absolute astrometric 
position. The size of G2's Br-7 emission was estimated by comparing the FWHM of G2 to that 
of nearby stars. To track the motions of G2 in the plane of the sky, its OSIRIS position was 
transformed into an absolute coordinate system. To construct an absolute reference frame, we 



used the pre cise locations of wel l characterized stars in the field (see figure 2(a) ), as measured with 
StarFinder (JDiolaiti et al.ll2000l ). The position of G2 from the two-dimensional Gaussian fit was 
added to the list of positions for each epoch. Then, these positions were matched to those from 
LGS AO K' observations of the same region at a nearby epoch using a second-order polynomial 
transformation to account for translation, rotation, and pixel scale differences between the imageqj. 
Once the positions in the OSIRIS frames were transformed into a K' epoch, they were transformed 
in to an absolut e refer ence frame in which SgrA*-radio is at rest at the origin, as originally described 
in lYelda et al.l (120101 ) and updated in Yelda et al. in prep. We investigated possible systematic 
astrometric effects in OSIRIS by cross-checking the OSIRIS astrometry of SO-2 to our standard 
NIRC2 analysis. The results are consistent, validating the use of OSIRIS astrometry. 



3.2. NIRC2 Imaging Measurements 

The L' position of G2 was estimated from calibrated near-infrared images in a similar fashion 
to that described above for the Br-7 nne maps. On e difference was that we deconvolved the 
indiv idual images using the Lucy- Richardson algorithm (JRichardsonlll972l : lLucylll974l : iLucy Sz Hook 
19920 in order to help isolate the point sources from the extended L' dust emission. Both the 
large-scale background and the PSF were estimated using StarFinder. After beam-restoring the 
individual frames with a Gaussian having FWHM of half the theoretical resolution limit (40 mas) , 
the frames were averaged to create a final image for each epoch (see figure |2(d)[ ) . The results of 
this analysis do not vary if the image is restored using a wider Gaussian beam. To locate stars used 
as astrometric reference sources, StarFinder was run on the combined (deconvolved) frame using 
the beam-restoring Gaussian as the reference PSF. G2's position was determined from an elliptical 
two-dimensional Gaussian fio The fitting of an elliptical Gaussian allows for extended structure 
and is analogous to the method used for the Br-7 detection. The positions for G2 from this method 
are consistent with those produced via StarFinder. The size of G2 was measured by comparing to 
the PSF of stars with a similar flux. The observed L' positions were transformed into an absolute 
coordinate system as described for the Br-7 detections. 



In 2012, only a linear transformation was used because fewer reference stars were available due to the smaller 
field of view (7 versus 19-27 in other epochs). 

2 In 2012, a two-component elliptical Gaussian was fit since G2 is blended with another point source. 



One complication for measuring the l! position of G2 in the plane of the sky is that G2 
was superimposed on a filament of dust emission during the early- to mid- 2000's, which caused 
significant astrometric bias (see figure [2(d)] ) , To quantify this effect, we performed a series of star- 
planting simulations using a 2012 May L' image, in which G2 is well off the filament. In each of the 
four simulations, an extra point source [V = 14) was planted at the expected G2 position at the 
time of earlier L' observations (2003.440, 2004.567, 2005.580, and 2006.385) based upon the orbit 
fit to the Br-7 measurements. The background dust emission does not have noticeable motions 
in the plane of the sky, thus the astrometric bias at these locations in the 2012 image should be 
representative of the bias at earlier epochs. The position of the artificial point source was extracted 
using StarFinder. 

To constrain the s pectral energy distr ibution of the obj ect, we measu r ed th e integrated flux of 
G2 i n the L f frame (e . g.lWitzel et al.ll2012l ). In contrast with lEckart et al.1 (J2013J), but in agreement 
with iGillessen et al.1 ( 20121 ). we detect no K' counterpart for G2; therefore, we performed star 
planting simulations to determine the K' upper limit. We used the Br-7 astrometry of G2 to 
predict the positions of G2 in the K' data. Then, an artificial point source was planted at G2's 



predicted position in the May 2010 image, when G2 is most isolated (see figure 2(c)). Using a 
modified version of the StarFinder algorithm (Boehle et al. in prep), we searched for a point 
source within a 3 pixel box centered on the planted star, given the PSF originally extracted from 
the May 2010 image. 



3.3. Orbital Modeling 

We fit an orbit to G2' s astrometric po s itions and R V measurements. L'-astrometric and RV 
measurements reported by IGillessen et al.l (120121 . 120131 ) were also included; however, we fit the 
Br-7 an d L' astrometry separately. We assumed G2 follows a purely Keplerian orbit around the 
central black hole. 

The orbit inference problem aims to find the most likely orbital parameters (eccentricity, period, 
time of periapse passage, inclination, position angle of the ascending node, and the longitude of 
periapse) given the data. In addition to G2's orbital parameters, the parameters of the gravitational 
potential need to be inferred (mass and distance to black hole, and its perceived 2D position and 
3D velocity). Since the observed motion of G2 does not entail enough information to constrain 
all these parameters, the orbit of SO-2, which has under gone a full orbit, was fit simultaneously 
and effectively determines the gravitational potential (see lMever et al.ll2012l ). Compared to earlier 
publications, we amend ed our Keplerian orbit fitting code to use the Bayesian sampler MultiNest 
(JFeroz &: Hobsonl 120081) . We y erified the results are equivalent with our previous Monte Carlo 
approach (e.g. iGhez et al.ll2008l ). 
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Fig. 1. — The evolving spectra of G2. The highly redshifted Br-7 emission line has a FWHM which 
increases with time (see table [1]). The Br-7 l me peak has been scaled such that it is the same height 
for each spectrum. Additionally, the spectra have been smoothed using a boxcar average with a 
width of 3 OSIRIS channels. 



4. Results 

Six key results emerge from our analysis of G2: 

» Apparent positional offset between the Br-j line emission and L' continuum emission of G2: 
The V astrometry is not consistent with our Br-7 astrometry, which is reported here for the 
first time (see figure [3]). The positional difference can be explained by astrometric bias in the 
V astrometry. In our star planting analysis, we find the extracted position of an artificial 
Jj source on the dust filament where G2 is located at early epochs differs from the input 
positions by an average of ~1.2 pix (12 mas). This is almost an order of magnitude larger 
than typical centroiding uncertainties at L' and the astrometric bias found in a region isolated 
from extended background flux. We therefore exclude L' astrometry from our primary orbital 
fits. 

» Revised Orbital Parameters: The Kepleri an orbit derived usin g Br-7 astrometry differs from 
the orbit derived using V astrometry (e.g. lGillessen et al.ll2013l ) by more than can be explained 



by formal measurement uncertainties (see figure [3]). We list the most likely orbital parameters 
from Br-7 i n table [2J This new orbit for G2 has pushed the closest approach date to March 
2014, has a closer periastron (130 AU, 1900^) and a higher eccentricity (0.9814 ± 0.0060). 

G2's compact size: G2's observed size is comparable to the Keck angular resolution in al- 
most all our observations; however, it exhibits marginal spatial extent in our highest quality 
spectroscopic data set. At this time (2006), the observed half-width at half-maximum is 
20 ± 4 mas after accounting for the observed PSF of SO-2 (FWHM = 74 ± 3 mas). G2 also 
shows a velocity gradient along its orbit, with the most highly redshifted portion of the object 
closest to the black hole, consistent with what is expected f or a tidal intera c tion e vent. The 



observed spatial extent is consistent with that presented in iGillessen et al.l (120121 ) , thus the 
Br-7 spatial extent of G2 along the direction of motion is relatively constant and corresponds 
to a size on the order of ~100 AU. 

Spectral Evolution: We find an increasing Br-7 nne width for G2, with FWHM measurements 
ranging from 135 ± 30 km s" 1 in 2006 to 706 ± 250 km s" 1 in 2012 (see table CD). These 



measurements are consistent with those presented by IGillessen et al.l ([2012J, |2013j) . 



Constant Brightness: Both at V and Br-7 the flux of G2 is constant within uncertainties. 
The integrated flux of G2 at V has an average value of 14.9 ± 0.3 mag. The Br-7 nne fluxes 
in 2006, 2009, 2010, and 2012 are 0.884 ±0.116%, 0.919 ±0.098%, 1.13 ±0.41%, 1.14 ±0.87% 
of S0-2's continuum flux between 2.17-2.18 //m, respectively. 

Source Color: Our star planting simulations indicate that the upper K' magnitude limit of 
G2 is K' = 20 and is ~ 2 magnitudes deeper than the earlier upper K' limit imposed on G2 



(JGillessen et al.ll2012l ). The deeper K' limit yields a dust temperature below ~500K if G2 is 



a pure gas cloud. If G2 has an underlying stellar source, a K' magnitude limit places distinct 



limits on the luminosity of the star. This limit is si milar to the expected K' magnitude for a 
low-mass T Tauri star in the Galactic Center (e.g. IScoville &: Burkertll2013l ). and eliminates 
a more massive star unless it is shrouded by dust which self-extincts the star at K' . 



5. Discussion 



The two orbits derived by iGillessen et al.l (|2012l . 120131 ) differ from each other by more than can 
be explained with formal measurement uncertainties. Systematic errors in the L' astrometry natu- 
rally explain the discrepancy. Because the offset becomes less significant and the RV measurements 
play an increasingly dominant role in the orbital modeling, the measurements seem to converge on 
the solution obtained with the less biased Br-7 astrometry and RV measurements. 

There are several implications of the new orbital model parameters. First, the revised periapse 
date of March 2014 and the closer periastron are important for the design and interpretation of 
SgrA* monitoring programs designed to test for an increased accretion flow. With a special NRAO 
call for proposals and tremendous attention called to this object, which was originally anticipated 
to experience periapse in June 2013, more than 30 programs have been approved for the summer of 
2013 covering radio to 7-ray wavelengths. While some models predict that a cloud would generate 



enhanced radio emissi on from SgrA* well in advance o f periapse passage ( e.g. iNarayan et al. 



2012 



Sadowski et al.ll2013al )bl). none has yet been detected (IKassim et al.ll2013l ). although this may not 
yet be expected with the revised periapse passage. 

Second, the geometry of the new orbit may call two previously suggested associations into 
question. Most importantly, the previous claimed "tail" of low surface brightness emission, which 
was central to the original claim that G2 is a pure gas cloud, is no longer securely affiliated with 
the compact "head", which has been analyzed in this work. The new orbit also falls 3.5a off 
the plane of the disk of young stars orbiting the central black hole at ~0.05 - 2 pc (JGenzel et al. 



Table 2. Orbital elements for G2 a 



Parameter 



Value b 



Time of closest approach (To ) 

Eccentricity (e) 

Period (P) 

Angle to periapse (ui) 

Inclination (i) 

Position angle of the ascending node (O) 



2014.21 ± 0.14 yrs 

0.9814 ± 0.0060 

276 ±111 yrs 

88 ± 6 dog 

121 ± 3 deg 

56 ± 11 deg 



"•Parameters des cribing the gravitational potential are found in 
iMever et all feoiizft . 

b Values provided are the mean and standard deviation of the 
marginalized one-dimensional distributions. 
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(a) A collapsed OSIRIS cube from May 2006 showing 
the continuum sources though the narrow-band Kn3 
filter. The position of G2, which has no detectable K 
continuum, is indicated with an arrow. 




(b) The Kn3 image shown in (a) is continuum- 
subtracted such that only Br-7 emission at G2's red- 
shift is shown. G2 is isolated in this continuum- 
subtracted map, minimizing the effect of stellar con- 
fusion on positional measurements. 




(c) A NIRC2 K' image from May 2010 with the 
OSIRIS position of G2 indicated by a 3<r contour. For 
reference, the detections of other K' sources are shown. 
This image was used to derive G2's K' magnitude up- 
per limit since it is most isolated in this imaging epoch. 




(d) A NIRC2 L' image from July 2012 which has been 
deconvolved and restored with a Gaussian Beam. The 
background estimation from StarFinder has also been 
added back in to show the dust filament which biases 
G2's L' position in early epochs. This feature extends 
upwards from the lower left of the image. G2 is blended 
with another point source in 2012. 



Fig. 2. — Images of the Galactic Center showing the Br-7 an d L' detections of G2, as well as the 
map in which the K' upper limit was derived. 
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20001 : iLevin k Beloborodovl 12003 : IPaumard et al.1 12009 : iLu et al.1 12003 : Bartko et al.1 12009 : Yelda 
et al. in prep) which has bee n suggested as the possible origin of the gas cloud through colliding 
stell ar winds of young stars (lOzernov et al.l 1 19971 : ICuadra et al.l 120061 ) or a central stellar source 
(e.g. iMurrav-Clav k Loeblboid ). 



Third, the even higher eccentricity provides strong constraints on G2's origin. It is notable that 
a high eccentricity is consistent with the outcome of three body exchanges between a binary system 
and the central black hole which may explain the dense concent ration of B stars (the S-stars) in the 



innermost regions of the Galaxy (e.g. lAlexander k Livid 120041 ) . If a triple star system is invoked 



rather than a binary, G2 could possibly be the result of a recent merger between two components 
to produce the observed physical properties (red, compact, and at times marginally resolved) and 
an exchange with the central black hole with the remaining component to produce the observed 
orbital parameters (highly eccentric, short period). While we favor a stellar model, this question 
will be soon addressed, as G2 should remain intact through its periapse if this is indeed correct. 

The authors would like to acknowledge the invaluable feedback from conference attendees at 
workshops where this work has been presented (AAS 221st Meeting, Keck 20th Anniversary Cele- 
bration, Galactic Nuclei Ringberg Workshop). Support for this work was provided by NSF grant 
AST-0909218. Data presented herein were taken at the W. M. Keck Observatory. The W. M. 
Keck Observatory is operated as a scientific partnership among the California Institute of Tech- 
nology, the University of California, and the National Aeronautics and Space Administration. The 
Observatory was made possible by the generous financial support of the W. M. Keck Foundation. 
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Fig. 3. — The orbital fit using the Br-7 astrometry is shown in blue and the fit derived using 
the L' astrometry is shown in red. Both solutions use radial velocity fits determined from the 
Br-7 emission line. Astrometric points obtained from data presente d in this paper are indicate d 
with filled points, and open circles indicate measurements presented in lGillessen et al.l (|2012l . 120131 ). 
Both V data sets appear offset from the Br-7 astrometry. 



